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Dark stars powered by dark matter annihilation have been proposed as the first luminous sources 
in the universe. These stars are believed to form in the central dark matter cusp of low-mass mini- 
halos. Recent calculations indicate stellar masses up to ~ 1000 and/or have very long lifetimes. 
The UV photons from these objects could therefore contribute significantly to cosmic reionization. 
Here we show that such dark star models would require a somewhat artificial reionization history, 
based on a double-reionization phase and a late star-burst near redshift z ~ 6, in order to fulfill 
the WMAP constraint on the optical depth as well as the Gunn-Peterson constraint at z ^ 6. This 
suggests that, if dark stars were common in the early universe, then models are preferred which 
predict a number of UV photons similar to conventional Pop. Ill stars. This excludes 800 Mq 
dark stars that enter a main-sequence phase and other models that lead to a strong increase in the 
number of UV photons. 

We also derive constraints for massive as well as light dark matter candidates from the observed 
X-ray, gamma-ray and neutrino background, considering dark matter profiles which have been steep- 
ened during the formation of dark stars. This increases the clumping factor at high redshift and 
gives rise to a higher dark matter annihilation rate in the early universe. We furthermore estimate 
the potential contribution from the annihilation products in the remnants of dark stars, which may 
provide a promising path to constrain such models further, but which is currently still uncertain. 



I. INTRODUCTION 

Growing astrophysical evidence suggests that dark 
matter in the universe is self-annihilating. X-ray obser- 
vations from the center of our Galaxy find bright 511 keV 
emission which cannot be attributed to single sources 
P, [1], but can be well-described assuming dark matter 
annihilation Q. Further observations indicate also an 
excess of GeV photons Q , of microwave photons [s] , and 
of positrons @ . A common feature of these observations 
is that the emission seems isotropic and not correlated 
to the Galactic disk. However, there is usually some dis- 
crepancy between the model predictions and the amount 
of observed radiation, which may be due to uncertainties 
in the dark matter distribution, astrophysical processes 
and uncertainties in the model for dark matter annihila- 
tion 0. 

It is well-known that weakly-interacting massive dark 
matter particles may provide a natural explanation of the 
observed dark matter abundance [1, Q . Calculations by 
Ahn et al. W\ indicated that the extragalactic gamma- 
ray background cannot be explained from astrophysical 
sources alone, but that also a contribution from dark 
matter annihilation is needed at energies between 1-20 
GeV. It is currently unclear whether this is in fact the 
case or if a sufficient amount of non-thermal electrons in 
active galactic nuclei (AGN) is available to explain this 
background radiation [ll|. Future observations with the 
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FERMI satellite [HI will shed more light on such ques- 
tions and may even distinguish between such scenarios 
due to specific si gna tures in the anisotropic distribution 
of this radiation [13 |. 

The first stars have been suggested to have high masses 
of the order ^ 100 Mq, thus pro viding powerful ionizing 
sources in the early universe [l3l [l3 | . The effect of dark 
matter annihilation on the first stars has been explored 
recently in different studies. Spolyar et al. [lB| showed 
that an equilibrium between cooling and energy deposi- 
tion from dark matter annihilation can always be found 
during the collapse of the proto-stellar cloud. This has 
been explored further by locco [l^ and Freese et al. [l3] , 
who considered the effect of scattering between baryons 
and dark matter particles, increasing the dark matter 
abundance in the star. locco et al. |18| considered dark 
star masses in the range 5 < < 600 Mq and cal- 
culated the evolution of the pre-main-sequence phase, 
finding that the dark star phase where the energy in- 
put from dark matter annihilation dominates may last 
up to lO'' yr. Freese et al. [l^ examined the formation 
process of the star in more detail, considering polytropic 
equilibria and additional mass accretion until the total 
Jeans mass of ~ 800 Mq is reached. They find that this 
process lasts for 5 x 10^ yr. They suggest that dark 
stars are even more massive than what is typically as- 
sumed for the first stars, and may be the progenitors for 
the first supermassive black holes at high redshift. locco 
et al. [lH, Taoso et al. [l^l and Yoon et al. 21 1 have 
calculated the stellar evolution for the case in which the 
dark matter density inside the star is enhanced by the 
capture of addition WIMPs via off-scattering from stel- 
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lar baryons. locco et al. fl§\ followed the stellar evolution 
until the end of He burning, Yoon et al. [2l[ until the end 
of oxygen burning and Taoso et al. [l^ until the end of H 
burning. Yoon et al. 21] also took the effects of rotation 
into account. The calculations found a potentially very 
long lifetime of dark stars and correspondingly a strong 
increase in the number of UV photons that may con- 
tribute to reionization. Dark stars in the Galactic center 
have been discussed by Scott et al. [2I, [2^ . 

Such models for the stellar population in the early 
universe imply that the first luminous sources produce 
much more ionizing photons, and reionization starts ear- 
lier than for a population of conventional Pop. HI stars. 
In fact, we recently demonstrated that reionization based 
on massive Pop. HI can well reproduce the observed 
reionization optical depth [23| . Increasing the number of 
ionizing photons per stellar baryon may thus reionize the 
universe too early and produce a too large reionization 
optical depth. This can only be avoided by introducing 
a transition to a stellar population which produces less 
ionizing photons, such that the universe can recombine 
after the first reionization phase. We therefore consider a 
double-reionization scenario in order to re-obtain the re- 
quired optical depth. We discuss such models in i ^IIII and 
demonstrate that some models of dark stars require con- 
siderable fine-tuning in reionization models in order to 
be compatible with the reionization optical depth from 
the WMAP 5-year data [p, HI] and to complete 

reionization at redshift z ~ 6 [27|. In mVl we show how 
such scenarios can be tested via 21 cm measurements. 

A further consequence of the formation of dark stars is 
the steepening of the density profiles in minihalos [16, 19], 
thus increasing the dark matter clumping factor with re- 
spect to standard NEW models. In ^Vl we estimate the 
increase in the clumping factor during the formation of 
dark stars and compare the calculation with our expecta- 
tion for conventional NEW profiles and heavy dark mat- 
ter candidates. In WI[ we perform similar calculations 
for the light dark matter scenario. Further discussion 
and outlook is provided in Will 



II. THE MODELS 

As discussed in the introduction, various models have 
been suggested for dark stars. The main difference be- 
tween these models comes from considering or neglecting 
scattering between dark matter particles and baryons. 
In addition, it is not fully clear how important a phase 
of dark matter capture via off-scattering from baryons 
actually is, depending on further assumptions on the 
dark matter reservoir. In the following, we will thus dis- 
tinguish between main-sequence dominated models and 
capture-dominated models. 



A. Main-sequence dominated models 

After an initial phase of equilibrium between cooling 
and heating from dark matter annihilation [Tsl . [T^ . [T^ , 
the dark star will contract further while the dark mat- 
ter annihilates away and the heating rate thus decreases. 
This duration of this adiabatic contraction (AC) phase is 
currently controversial: While locco et al. pJi] find it to 
be in the range of (2 — 20) x 10^ yr, Freese et al. [l^ re- 
quire about 10® yr. However, with a surface temperature 
of ^ 6000 K, the stars are rather cold in this phase, and 
thus will not contribute significantly to reionization. The 
uncertainty in the duration of the AC phase is therefore 
not crucial in this context. 

If the elastic scattering cross section as well as the dark 
matter density around the star are sufficiently large, the 
star will enter a phase which is dominated by the capture 
of further dark matter particles. Such a scenario will be 
discussed in more detail in the next subsection. Here, we 
assume that the elastic scattering cross section is either 
too small, or that the dark matter reservoir near the star 
is not sufficient to maintain the capture phase for long. 
Then, the star will enter the main-sequence phase (MS), 
in which the luminosity is generated by nuclear burning. 
Stars with ~ 1000 Mq are very bright in this phase, 
and emit ^ 4 x 10* hydrogen- ionizing photons per stellar 
baryon during their lifetime [1^ [2^. We will refer to 
stars of such type, which have only a short or even no 
phase driven by dark matter capture, as MS-dominated 
models. 

Eor the case of MS-dominated models, we will focus 
essentially on the very massive stars suggested by Ereese 
et al. [1^. For stars in the typical Pop. Ill mass range, it 
has been shown elsewhere [e. g. [IJ] that they are consis- 
tent with reionization constraints. A star with ^ 800 Mq 
forming in a dark matter halo of ~ 10® Mq corresponds 
to a star formation efficiency of 1%, which we adopt for 
this case. 



B. Capture-dominated models 

Eor a non-zero spin-dependent scattering cross section 
between baryons and dark matter particles, stars can 
capture additional WIMPs which may increase the dark 
matter density inside the star. Eor a cross section of the 
order 5 x 10^"^^ cm^ and an environmental dark mat- 
ter density of ^ 10^" GeV cm^'^, this contribution be- 
comes significant and alters the stellar evolution during 
the main sequence phase. We will refer to such a scenario 
as a capture-dominated (CD) model. These phases have 
been studied in detail by locco et al. [l^, Taoso et al. 
[20| and Yoon et al. [2l|. They found that the number of 
ionizing photons produced by such stars may be consid- 
erably increased with respect to high-mass stars without 
dark matter annihilation effects, which is mostly due to a 
longer lifetime. In particular for dark matter densities of 
(1 — 5) X 10^" GeV cm~'^ , the number of produced ioniz- 
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ing photons may be increased by up to two orders mag- 
nitude, while it decreases rapidly for larger dark matter 
densities, and the number of ionizing photons per baryon 
even drops below the value for Pop. II stars at. threshold 
densities of 1 x 10^^ GeV cm^^. As Yoon et al. 2l'] found 
only a weak dependence on stellar rotation, we will not 
explicitly distinguish between models with and without 
rotation in the follow. 

For the calculation of reionization, we will focus on 
some representative models of Yoon et al. [2l| in the fol- 
lowing. However, we point out that there are still signifi- 
cant uncertainties in these models, in particular the dark 
matter parameters and the lifetimes of the stars. The lat- 
ter should be seen as upper limits, as they assume that 
a sufficient reservoir of dark matter is available in the 
stellar neighborhood to allow for ongoing dark matter 
capture. This may however be disrupted by dynamical 
processes. An apparent disagreement of dark stars in the 
early universe with our reionization model may thus in- 
dicate that the stellar lifetimes are indeed smaller due to 
such processes. 



III. REIONIZATION CONSTRAINTS 

In this section, we briefly review our reionization model 
and discuss reionization histories for main-sequence and 
capture-dominated models. These calculations implicitly 
assume annihilation cross sections of the order 10~^^ cm^ 
and dark matter particle masses of the order 100 GeV, 
the values which are typically adopted in dark star mod- 
els. In such models, dark matter annihilation does not 
contribute to cosmic reionization [23 |. The chemistry 
in the pre-ionization era is thus unchanged and well- 
described by previous works [13, [Ml, [H, [sl] , such that 
the initial conditions for star formation are unchanged. 
Considering higher annihilation cross sections essentially 
yields an additional contribution to the reionization op- 
tical depth, which would sharpen the constraints given 
below. 



A. General approach 

Our calculation of reionization is based on the frame- 
work developed by Schleicher et al. [241. w hich we have 
implemented in the RECFAST code [ill [H, [3|. We 
will review here only those ingredients which are most 
relevant for this work. During reionization, the IGM con- 
sists of a two-phase medium, i. e. a hot ionized phase and 
a rather cold and overally neutral phase. The relative size 
of these phases is determined from the volume-filling fac- 
tor Qh+ of the H+ regions [li S E E3 as a 
function of redshift, given by 



dz 



H{z){l + z) 



(1) 



where C{z) = 27 Am exp(-0.114z + 0.001328z2) is the 
clumping factor [i^ , ne,H+ the number density of ionized 
hydrogen, a a the case A recombination coefficient j43l |. 
H{z) the Hubble function, tih the mean neutral hydro- 
gen density in regions unaffected by UV feedback and 
duph/dz the UV photon production rate. Our model 
consists of ordinary differential equations (ODEs) for the 
evolution of temperature T and ionized fraction Xi in the 
overall neutral medium. For the application considered 
here, the dominant contribution to the effective ionized 
fraction x^f / = Qh+ + (1 ~ Qn+ja^i and the effective tem- 
perature Teff — 10^ K (5h+ + T{1 — Qh+) comes indeed 
from the UV feedback of the stellar population, i. e. from 
the hot ionized phase. According to Gnedin and Hui [i^l 
and Gnedin [45|, we introduce the filtering mass scale as 



M 



2/3 



1/2- 



(2) 



where a = (1 + z) ^ is the scale factor and Adj the ther- 
mal Jeans mass, given as 



Mj = 2Me 
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Here, Cs is the sound speed evaluated at temperature 
Tgff, in order to take into account the backreaction of 
heating on structure formation. In this framework, the 
production of UV photons can be described as 



dnph/dz ^ ^ rf/coii 
rtjj dz ' 



(4) 



where e = AnoA/cscMon, with Aho = 4./{A-3Yp) = 1.22, 
A'ion the number of ionizing photons per stellar baryon, 
/* is the star formation efficiency and /esc the escape frac- 
tion of UV photons from their host galaxies. The quan- 
tity /cou denotes the fraction of dark matter collapsed 
into halos, and is given as 



/coll = erfc 



%/2cr(M„iin) 



(5) 



where Mmi„ = min(Mi;^, 10^ M©), 5c = 1.69/D{z) is the 
linearized density threshold for collapse in the spherical 
top-hat model and <7{AIyain) describes the power associ- 
ated with the mass scale Mmin- 

A relevant question in this context is also the role of 
Lyman- Werner (LW) feedback, which may suppress the 
star formation rate in low-mass halos. The role of such 
feedback has been addressed using different approaches. 
For instance, Machacek et al. [46|, O'Shea and Norman 
[13] and Wise and Abel (isj have addressed this ques- 
tion employing numerical simulations in a cosmological 
context, assuming a constant LW-background radiation 
field. These simulations indicated that such feedback can 
delayed star formation considerably. 

More self-consistent simulations show, however, that 
the above calculations overestimated the role of LW- 
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FIG. 1: The evolution of the effective ionized fraction aj^g, for 
reionization models with main-sequence dominated dark stars 
(see Table IB). Models MS 1 and MS 2 can be ruled out by 
reionization constraints, while models MS 3 and MS 4 require 
a sudden increase in the star formation rate by a factor of 
30 at redshift 6.5. It appears more realistic to assume lower 
masses and star formation efficiencies to reconcile dark star 
models with observations. 



feedback. Considering single stellar sources and neglect- 
ing self-shielding, Wise and Abel [4§| showed that LW- 
feedback only marginally delays star formation in halos 
that already started collapsing before the nearby star 
ignites. More detailed simulations taking into account 
self-shielding show that the star formation rate may be 
changed by only 20% in the presence of such feedback 
[50| . This is due to the rapid re-formation of molecular 
hydrogen in relic HII regions, which leads to abundances 
of the order 10^'*. Such abundances effectively shield 
against LW- feedback and make it ineffective [5l|. This 
is the point of view adopted here, which may translate 
into an uncertainty of ~ 20% in the star formation rate. 
In fact, in scenarios involving dark matter annihilation, 
H2 formation and self-shielding could be even further en- 
hanced compared to the standard case [52 1. 

The models have to reproduce the reionization optical 
depth given by r = 0.087± 0.017 [2^ and fully ionization 
at z ^ 6 27]. In the following, we will try to construct 
appropriate reionization histories for the different dark 
star models. 



B. Reionization with MS-dominated dark stars 

As shown previously [l^], MS-dominated dark stars 
with ~ 1000 Mq would significantly overproduce the 
reionization optical depth if this type of stars had been 
common throughout the early universe. If, on the other 
hand, MS-dominated stars only had mass scales of ~ 
100 Af0, comparable to conventional Pop. Ill stars, reion- 
ization could not discriminate between them and conven- 
tional Pop. Ill stars, and dark stars would be compatible 



with observations. Alternatively, as explained in the in- 
troduction, a transition in the stellar population might 
help to alleviate the problem for high-mass dark stars. 
We will explore this possibility in more detail to work 
out whether such a scenario is conceiveable. 

Numerical simulations by Dove et al. [s^ , Ciardi et al. 
[5^ and Fujita et al. 55] indicated rather high es- 
cape fractions of order 100% for massive Pop. Ill stars. 
Wood and Loeb [HH] found rather low escape fractions 
below 10%, while radiation hydrodynamics simulations 
by Whalen et al. [13] show that such stars can eas- 
ily photo-evaporate the minihalo. Here we adopt the 
point of view that indeed massive stars can photoe- 
vaporate small minihalos, but that the escape fraction 
will be reduced to ~ 10% in atomic cooling halos that 
have virial temperatures larger than 10"* K. Thus, we set 

/esc = 1 if the filtering mass is below the mass scale 

/ \ 3/2 

Mc = 5 X lO^M© ( j that corresponds to the virial 

temperature of 10"^ K [5l,[5§], and /esc = 0.1 in the other 
case. To reflect the expected stellar mass of ~ 800 Mq, 
we choose a star formation efficiency of ^ 1%, an 
order of magnitude higher than what we expect for con- 
ventional Pop. Ill stars |.24i] . 

Assuming that reionization is completely due to these 
MS-dominated dark stars (model MS 1), we find that the 
universe is fully ionized at redshift Zrcion = 15.5 and the 
reionization optical depth is Tj-oion ~ 0.22, i. e. signifi- 
cantly larger than the WMAP 5 optical depth (see Fig. 
[T]). Such a model is clearly ruled out. 

To reconcile the presence of such massive dark stars 
with observations, one could invoke a double-reionization 
scenario, assuming a transition to a different mode of star 
formation induced by the strong UV feedback of MS- 
dominated dark stars. In fact, even for conventional star 
formation models, it is discussed that such UV feedback 
may lead to a less massive mode of star formation (goI , l6ll , 
[g^ . In addition, chemical enrichment should facilitate 
such a transition as well 41, 63, 64, 65, 66, 67], although 
it is unclear how well metals will mix with the pristine 
gas. We assume that the transition to a low-mass star 
formation mode with a Scalo-type IMF [g^] happens at 
redshift 15.5, when the universe is fully ionized and UV 
feedback fully effective. For the subsequent Pop. II stars, 
we assume a star formation efficiency of = 5 x 10^'^ 
and A^ion = 4 x 10'^ UV photons per stellar baryon. 

Corresponding photon escape fractions are highly 
uncertain. Observations of Steidel et al. [i^ indicate 
an escape fraction of 10% at z ~ 3, while others 
find detections or upper limits in the range 5 — 10% 
0, [Zil, Izi Izi- We adopt the generic value of 10% for 
simplicity, though our results do not strongly depend on 
this assumption. For this scenario, to which we refer as 
model MS 2, we find an optical depth Troion = 0.082 well 
within the WMAP constraint, but the universe does not 
get fully ionized until redshift zero. This scenario is thus 
rejected based on the constraint from quasar absorption 
spectra [271]. 
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To fulfill both the WMAP constraint as well as 
full- ionization at z '--^ 6, we need to introduce an addi- 
tional transition in our model. At redshift Zburst — 6.5, 
we increase the star formation efficiency to 15%. This 
might be considered as a sudden star burst and results 
in fuU-ionization at z = 6.2. In this case, we find 
Trcion — 0.116, which is within the 2a range of the 
WMAP data. However, we are not aware of astrophysi- 
cal models that provide a motivation for such a sudden 
star burst that increases the star formation rate by a 
factor of 30. Based on gamma-ray burst studies, Yiiksel 
et al. [zl] showed that the cosmic star formation rate 
does not change abruptly in the redshift range between 
redshift zero and Zburst — 6.5. Such a sudden burst is 
thus at the edge of violating observation constraints. 

To improve the agreement with WMAP, one can con- 
sider to shift the first transition to zpop ii — 18 where 
full ionization is not yet reached (model MC 4), which 
yields an optical depth Tioion = 0.086, in good agreement 
with WMAP. At this redshift, 68% of the universe are 
already ionized, so UV feedback might already be active 
and induce a transition in the stellar population. The 
results are given in Fig. [1] and summarized in Table [H 



Model 


^Pop II 


■^burst 


^icion 




MS 1 






0.22 


15.5 


MS 2 


15.8 




0.078 


never 


MS 3 


15.5 


6.5 


0.116 


6.2 


MS 4 


18. 


6.5 


0.086 


6.2 



TABLE I: Reionization models for MS-dominated dark stars. 
The parameters zpop ii and ^burst give the transition redshifts 
to a mode of Pop. II star formation and to the sudden star 
burst, while Trcion is the calculated reionization optical depth 
and Zf the redshift of full ionization. 

However, we find that only models MS 3 and MS 4 
cannot be ruled out observationally. These models re- 
quire two severe transitions in the stellar population and 
cannot be considered as "natural". Improved measure- 
ments of the reionization optical depth from Planck [l25j | 
will remove further uncertainties and may rule out model 
MS 3 as well. From a theoretical point of view, it must 
be checked whether strong UV feedback can lead to the 
required transition to a low-mass star population, and in 
addition, the plausibility of a sudden star burst near red- 
shift 6 must be examined as well. In summary, it seems 
more plausible to conclude that MS-dominated dark stars 
were less massive than suggested b y F reese et al. [T^, as 
already hinted by Schleicher et al. [24j . 



C. Reionization with CD darlc stars 

For CD dark star models, the situation is complicated 
by the fact that the number of UV photons per stellar 
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FIG. 2: The evolution of the effective ionized fraction x^^g, for 
reionization models with capture-dominated dark stars (see 
Tableinil. Models CD la, CD lb, CD 2a, CD 2b and CD 3 are 
ruled out due to reionization constraints, while the remaining 
models require an artificial star burst. 



baryon, iVion , is model-dependent and changes with the 
environmental dark matter density, px. We select three 
representative models of Yoon et al. [2l| , which assume a 
spin-dependent scattering cross section of 5 x 10~^^ cm^ 
(see Table HIl . In general, stellar models depend on the 
product of this scattering cross section with the threshold 
dark matter density at the stellar radius [13, [2l| . Lower 
elastic scattering cross sections therefore correspond to 
going to smaller threshold densities at the same elastic 
scattering cross section. 

In the models CD 1 and 2, iVjon is larger than for con- 
ventional Pop. HI stars, while in the model CD 3, it is 
even less than in the case of Scalo-type Pop. II stars. 
Such a low luminosity is unlikely to photo-evaporate star- 
forming halos, and we thus adopt /esc ~ 10% for this 
case. However, such Scalo-type Pop. II stars are ruled out 
as sole sources for reionization |24| . As we show in Fig. 
[21 even with a high star formation efhciency of = 1%, 
they never ionize the universe completely. 

In principle, one could consider the presence of other 
sources to ionize the universe. While dark stars of type 
CD 3 may be the first stars to form, one might envision a 
transition to a stellar population with the power to ion- 
ize the universe. This transition is unlikely due to UV 
feedback, as UV feedback from dark stars is rather weak 
in this scenario. One thus has to rely on effective mixing 
of the produced metallicity, or assume that the first stel- 
lar clusters in atomic cooling halos contain a sufficient 
number of massive stars to reionize the universe [63j . 

For the other two models, TVion is significantly larger 
and we adopt the procedure from the previous subsection, 
such that /osc depends on the filtering mass. We adopt a 
star formation efficiency of /* = 0.1%. We examine the 
reionization models given in Table |TT1 which essentially 
follow the philosophy of the models from the previous 
section. We calculate the reionization history for the case 



6 



where these dark stars are sole sources (CD la, CD2a) 
and find that the optical depth is considerably too high. 
We then determine the redshift where the universe is fully 
ionized and assume a transition to Pop. II stars at this 
redshift. In addition, to obtain full ionization at redshift 
6, we assume a late star burst as in the models MS 3 
and MS 4. This approach corresponds to the models CD 
lb and CD 2b, and yields optical depth that are at least 
within the 2a error of WMAP 5. In the models CD Ic 
and CD 2c, we improve the agreement with WMAP by 
introducing the Pop. II transition at an earlier redshift. 

The results are given in Fig. [21 Again, it turns out 
that somewhat artificial models are required to allow for 
an initial population of CD dark stars. The best way to 
reconcile these models with the constraints from reion- 
ization might be to focus on those models that predict 
a parameter A'ion which is closer to the Pop. Ill value 
of 4 X 10^. This may be possible, as the transition from 
the models CD 1 and 2 to CD 3 is likely continuous, 
and an appropriate range of parameters may exist to re- 
concile models with observations. This would require a 
px between 10" GeV cm"^ and 10^^ GeV cm^^. As 
mentioned earlier, the apparent violation of reionization 
constraints by some models depends also on the uncer- 
tainties in the stellar lifetime. If the dark matter reservoir 
near the star is destroyed earlier due to dynamical pro- 
cesses, the lifetime may be significantly reduced. Also, we 
stress that the conclusions depend on the adopted elastic 
scattering cross section and the dark matter density in 
the environment. The discussion here is limited to those 
models that have previously been worked out in detail. 



IV. PREDICTIONS FOR 21 CM 
OBSERVATIONS 

While some of the models suggested above essentially 
co-incide with standard reionization by mimicing the ef- 
fects of conventional Pop. Ill stars, others may have a 
very distinctive signature, as they consist of a double- 
reionizat ion p hase, and upc oming 21 cm telescopes like 
LOFAR [Hi] or SKA [l2|] can thus verify or rule out 
such suggestions. The calculation shown in Fig. [3] is 
based on the double reionization model MS 4, but clearly 
the models MS 3, CD lb, CD Ic, CD 2b and CD 2c yield 
similar results. In such a double-reionization scenario, 
the gas is heated to ^ 10^ K during the first reionization 
epoch. Assuming that the first reionization epoch ends 
at redshift zpop n, the gas temperature in the non- ionized 
medium will then evolve adiabatically as 



T - 10'' K 



l + z 



1 + ^^Pop II 



(6) 



In addition, the previous reionization phase will have es- 
tablished a radiation continuum between the Lyman a 
line and the Lyman limit, where the universe is opti- 
cally thin, apart from single resonances corresponding to 
the Lyman series. This radiation is now redshifted into 
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FIG. 3: 21 cm signatures of double-reionization scenarios 
(here MS 4 from Table |l]). Given is the evolution after the 
first reionization phase, when the H gas is heated from the 
previous ionization. Top: HI gas temperature, here identi- 
cal to the spin temperature. Middle: Expected mean 21 cm 
brightness fluctuation. Bottom: Frequency gradient of the 
mean 21 cm brightness fluctuation. 



the Lyman series and may couple the spin temperature 



^as temperature T via 
|76| . In fact, a small 



Tspin of atomic hydrogen to the 
the Wouthuysen-Field effect [tsI . 
amount of Lyman a radiation suffices to set Tgpin = T 
[ttI . ItsI . [z§|, which we assume here. Also, as the uni- 
verse is optically thin to this radiation background, even 
Pop. II sources will suffice to couple the spin temperure 
to the gas temperature. The mean 21 cm brightness tem- 
perature fluctuation is then given as 



5Ti, = 27a;H(l + <5) 



0.023 



0.15 \ + z' 
Vl,nh? 10 



1/2 
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Reion. model 


px/10 


12 




/* 


2Pop II 


'^rcion 


CD la 


0.01 GeV cm" 


-3 


1.75 X 10^ 


0.1% 


- 


0.162 


CD lb 


0.01 GeV cm" 


-3 


1.75 X lO"^ 


0.1% 


12.7 


0.109 


CD Ic 


0.01 GeV cm" 


-3 


1.75 X 10^ 


0.1% 


14.5 


0.089 


CD 2a 


0.05 GeV cm" 


-3 


2.4 X 10*' 


0.1% 




0.283 


CD 2b 


0.05 GeV cm" 


-3 


2.4 X 10® 


0.1% 


21.6 


0.106 


CD 2c 


0.05 GeV cm" 


-3 


2.4 X lO** 


0.1% 


23 


0.084 


CD 3 


1 GeV cm" 


-3 


1.1 X 10^ 


1% 




0.004 



TABLE II: Reionization models for CD dark stars stars. The number of ionizing photons was determined from the work of 
Yoon et al. [21I]. The parameters zpop 11 and Zburst give the transition redshifts to a mode of Pop. II star formation and to the 
sudden star burst, while Tj-oion is the calculated reionization optical depth and Z[ the redshift of full ionization. The calculation 
assumes a spin-dependent scattering cross section of 5 x 10^"*^ cm^. As stellar models depend on the product of this cross 
section with the threshold dark matter density, the effect of a lower scattering cross section is equivalent to a smaller threshold 
density. 



Ts-Tr 
Ts 



H{z)/{l + z) 

dv\\/ dr\\ 



mK, 



(7) 



where denotes the neutral hydrogen fraction, 5 the 
fractional overdensity, fi;,, $7^ the cosmological density 
parameters for baryons and total matter, h is related to 
the Hubble constant Hq via h = 7Jo/(100km/s/Mpc), Tr 
the radiation temperature and dv\\/dr\\ the gradient of 
the proper velocity along the line of sight, including the 
Hubble expansion. We further calculate the frequency 
gradient of the mean 21 cm brightness temperature fluc- 
tuation to show its characteristic frequency dependence. 
In Fig. [31 we show the evolution of the gas temperature, 
the mean 21 cm brightness fluctuation and its frequency 
gradient for model MS 4. 

As pointed out above, we expect similar results for 
other double-reionization models because of the char- 
acteristic adiabatic evolution of the gas and spin 
temperature. The decrease of the spin temperature 
with increasing redshift is a unique feature that is not 
present in other models that like dark matter decay 
[80| or ambipolar diff'usion heating from primordial 
magnetic fields fsi], [s^, HI] , which may also increase the 
temperature during and before reionization. 



V. COSMIC CONSTRAINTS ON MASSIVE 
DARK MATTER CANDIDATES 

In typical dark star models, it is assumed that massive 
dark matter candidates like neutralinos with masses of 
the order 100 GeV annihilate into gamma-rays, electron- 
positron pairs and neutrinos [H, [H, [TtI. [Tsl. [T9l[2ll. [83| . 
Similar to the constraint on high-redshift quasars from 
the X-ray background [85, 86, 87], the gamma-ray and 
neutrino backgrounds allow to constrain the model for 
and the amount of dark matter annihilation. As de- 
tailed predictions for the decay spectra are highly model- 
dependent, it is typically assumed that roughly 1/3 of the 
energy goes into each annihilation channel. Constraints 
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FIG. 4; The predicted gamma-ray background due to di- 
rect annihilation into gamma-rays in the presence of adia- 
batic contraction during the formation of dark stars, and the 
background measured by EGRET (squares) [92]. One finds 
two peaks in the annihilation background for a given parti- 
cle mass: One corresponding to annihilation at redshift zero, 
and one corresponding to the redshift where the enhancement 
from adiabatic contraction was strongest. 



on such scenarios are available from the Galactic cen- 
ter and the extragalactic gamma-ray and neutrino back- 
grouns [H, [H, [13, mi • In this section, we consider how 
such constraints are affected when the increase in the 
annihilation rate due to enhanced dark matter densities 
after the formation of dark stars is taken into account. 



A. Gamma-ray constraints 

We adopt the formalism of Mack et al. [9l[ who re- 
cently addressed the direct annihilation of massive dark 
matter particles into gamma-rays. The background in- 
tensity In is given from an integration along the line of 



8 



sight as 



c / dzP^([l + z\v, z) 



i?(z)(l + z)4 



(8) 



where P^(v,z) is the (proper) volume emissivity of 
gamma-ray photons, which is given as 



Pi, = abS ((1 + z)!^ - todm) keV 



'71 



(9) 

where (av) = 3 x 10 cm'^ s ^ denotes the thermaUy- 
averaged annihilation cross section, at = 1/3 is the 
adopted branching-ratio to gamma-rays and ttidm the 
mass of the dark matter particle in keV. C-y refers to 
the dark matter clumping factor. This clumping factor 
depends on the adopted dark matter profile and the as- 
sumptions regarding substructure in a halo [H, [H, HI, 
l96j . Here we use the clumping factor for a NFW dark 
matter profile 1971 which has been derived by Ahn and 
Komatsu [sl, l98l |. For z < 20, it is given in the absence 
of adiabatic contraction as a power-law of the form 



Cdm = CDM(0)(l + z)- 



(10) 



where Cdm(O) is the clumping factor at redshift zero 
and /3 determines the slope. For a NFW profile [IJ, 
Cdm(O) ~ 10^ and f3 ^ 1.8. The enhancement due to 
adiabatic contraction is taken into account by defining 



C-f — C'DM/cnhj 



(11) 



where the factor /enh describes the enhancement of the 
halo clumping factor due to adiabatic contraction (AC). 
We have estimated this effect based on the results of locco 
et al. [13, comparing a standard NFW profile with the 
enhanced profile that was created during dark star for- 
mation. We only compare them down to the radius of the 
dark star and find an enhancement of the order ~ 10"^. 
For the NFW case, the clumping factor would be essen- 
tially unchanged when including smaller radii as well, 
while the AC profile is significantly steeper and the con- 
tribution from inside would dominate the contribution 
to the halo clumping factor. However, as the annihila- 
tion products are trapped inside the star, it is natural to 
introduce an inner cut-off at the stellar radius. In addi- 
tion, we have to consider the range of halo masses and 
redshifts in which dark stars may form. We assume that 
the halo mass must be larger than the filtering mass to 
form dark stars. However, there is also an upper mass 
limit. Halos with masses above 



Mc = 5 X lO^Mp 



10 

l + z 



3/2 



correspond to virial temperatures of 10^ K [sl] and are 
highly turbulent [s^. It seems thus unlikely that stars 
will form on the very cusp of the dark matter distribution 
in such halos, and more complex structures may arise. 
We thus assume that dark stars form in the mass range 
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FIG. 5: The predicted neutrino background due to direct an- 
nihilation into neutrinos in the presence of adiabatic contrac- 
tion during the formation of dark stars, and the atmospheric 
neutrino background j99j] . One finds two peaks in the annihi- 
lation background for a given particle mass: One correspond- 
ing to annihilation at redshift zero, and one corresponding to 
the redshift where the enhancement from adiabatic contrac- 
tion was strongest. 



between Mp and Mc- Once Mc becomes larger than 
Mp, dark star formation must end naturally. In fact, it 
may even end before, as discussed in §1111 To obtain the 
highest possible effect, we assume that dark stars form 
as long as possible. We thus have 



/cnh =1 + 10 



~l3./coll(AfF) — fcon{Mc) 

fconiMp) 



(13) 



In Fig. m we compare the results with EGRET observa- 
tions of the gamma-ray background [o^l ■ In the absence 
of adiabatic contraction, the predicted background peaks 
at the contribution from redshift zero [oij. We find that 
the enhancement of annihilation due to adiabatic con- 
traction produces a second peak in the predicted back- 
ground which originates from higher redshifts. In this 
scenario, particle masses smaller than 30 GeV can thus 
be ruled out. 



B. Neutrino constraints 

The contribution to the cosmic neutrino flux can be 
obtained in analogy to Eq. ([8]) . As recent works [H, [93] , 
we adopt an annihilation spectrum of the form 



(12) P, = atS {{1 + z)iy - muu) 



keV 



keV (cr-u)n^MC, 



ncutnnoi 



(14) 



which is analogous to the spectrum for annihilation into 
gamma-rays. The branching ratio to neutrinos is as- 
sumed to be 1/3 as well, and the annihilation comes 
from the same dark matter distribution, thus yielding 



C^. The atmospheric neutrino background 
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matter left over in the final remnant can contribute). In 
this case, we have a proper volume emissivity 
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FIG. 6: The maximum gamma-ray background due to direct 
annihilation into gamma-rays in the remnants of dark stars, 
and the background measured by EGRET '92]. The actual 
contribution to the gamma-ray background is highly model- 
dependent (see discussion in the text). 



has been calculated from different experiment s with gen- 
erally good agreement [9i,[ioO,[Ioiillol[lQl. loccolTel 
adopted a similar atmospheric neutrino flux for compari- 
son with the expected neutrino flux from dark stars. We 
adopt here the data provided by Honda et al. [9§| and 
compare them to the predicted background in Fig.O The 
predicted background is always well below the observed 
background. 



C. Emission from dark star remnants 

In the previous subsections, we have included the en- 
hancement of the halo clumping factor down to the stel- 
lar radius, as by definition the annihilation products on 
smaller scales are trapped inside the star. At the end of 
their lifetime, these stars may explode and the baryon 
density in the center may be largely depleted. The dark 
matter density has certainly been significantly reduced 
due to annihilations during the lifetime of the star, but it 
may still be enhanced compared to the usual NFW case. 
A detailed calculation of this effect is strongly model- 
dependent. As we have seen above, the strongest con- 
straints are obtained for direct annihilation into gamma- 
rays, which is the case we pursue here in more detail. 

So far, we assumed that dark stars form in halos be- 
tween the filtering mass Mp and the mass corresponding 
to a virial temperature of 10* K, Mc- To obtain an upper 
limit, it is sufficient to assume that in all halos above Mp 
a dark star remnant will form at some point. Such an as- 
sumption clearly overestimates the total contribution at 
low redshift. When the dark star has formed, a fraction 
/core ~ 10~^ of the dark matter from the total halo is in 
the star [l^. For the upper limit, we assume that the 
total amount of dark matter in star will contribute to 
the X-ray background (in fact, however, only the dark 



P„ = S ((1 + z)iy - uium) keVasii /,./a 

keV 



dt 



(15) 



where riDM is the mean proper number density of dark 
matter particles, todm the particle mass in keV and 
df coll /dt can be evaluated from Eq. ([5]). The model- 
dependent factor fr determines which fraction of the dark 
matter in the star will be left in the remnant. We adopt 
/r = 1 to obtain an upper limit. The factor /a deter- 
mines the fraction of the remaining dark matter which 
actually annihilates, which we set to /a = 1 as well. As 
in WI A| asii = 1/4 is the fraction of electron-positron 
annihilations per one dark matter annihilation process, 
corresponding to annihilation via positronium formation. 
In Fig. 6l we compare the results with EGRET obser- 
vations 921 ■ We find that the maximum contribution is 
clearly above the observed background. 

Whether this maximum contribution can be reached, 
is however uncertain and the previous work in the litera- 
ture only allows one to make rather crude estimates. For 
instance locco et al. [ll| calculate the density profile for a 
fiducial 100 Mq protostar, finding that the density within 
the star roughly scales with outside a plateau at a 
radius r ~ 10"'^^ cm. At the stellar radius of ~ 10^"* cm, 
the dark matter density is still ^ 10^^ GeV cm^"^. The 
timescale to remove this dark matter enhancement by 
annihilation is ~ 100 Myr for 100 GeV neutralinos. We 
need to estimate which fraction of the dark matter inside 
the star will be left at the end of its life, where the gas 
density is expelled by a supernova explosion and the dark 
matter annihilation from this region may contribute to 
the cosmic gamma-ray background. 

Yoon et al. [2l[ adopted a timescale of 100 Myr, the 
typical merger timescale at these redshifts, as the maxi- 
mum lifetime for dark stars. Depending on the scattering 
cross section and the environmental density, the actual 
dark star lifetime may be considerably shorter. Indeed, 
as we showed in § IIII C[ it is difficult to reconcile lifetimes 
of ~ 100 Myr with appropriate reionization scenarios. It 
is therefore reasonable to assume shorter timescales. In 
such a case, a reasonable estimate is that 40% of the 
dark matter inside the star would be left at the end of 
its life. This would still be enhanced compared to the 
standard NFW profile. In this case, the parameter fr 
is ~ 40%, and fa may be of order 1, as the annihila- 
tion timescale is comparable to the Hubble time. We 
note that these numbers are highly uncertain, in particu- 
lar regarding the exact evolution of dark matter density 
during the lifetime of the star, the effect of a supernova 
explosion on the dark matter cusp as well as the conse- 
quences of minor mergers. 

There is however also a viable possibility that the 
dark matter distribution inside the star is significantly 
steeper than assumed above. In the case of dark matter 
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capture by ofF-scattering from baryons, the dark matter 
density inside the star follows a Gaussian shape and is 
hi g hly con centr ated in a small region of r ~ 2 x 10^ cm 
(To. [20. [2ll . Il04l |. The implications are not entirely clear. 
If capture of dark matter stops at the end of the life of 
the star, the density inside the star will annihilate away 
quickly, and no significant contribution may come from 
the remnant. If, on the other hand, dark matter cap- 
ture goes on until the end of the life of the star, a con- 
tribution to the background seems viable. In summary, 
this may provide a potential contribution to the cosmic 
background, but its strength is still highly uncertain and 
should be explored further by future work. 



D. Dependence of dark star models on the 
neutralino mass 

We conclude this section with a discussion on the con- 
straints from cosmic backgrounds for different neutralino 
masses. As dark star models in the literature mostly 
consider neutralinos of 100 GeV, there are uncertain- 
ties that need to be addressed when considering differ- 
ent neutralino masses. For models involving the capture 
of dark matter, locco (IS] states that the mass of the 
neutralino does not change the annihilation luminosity. 
Taoso et al. [2^ find that variations due to different neu- 
tralino masses are less than 5%. While these results may 
hold for high masses, Spergel and Press |l05j | showed that 
for neutralino masses below 4 GeV, they would evaporate 
from the star, as scattering with baryons can upscatter 
them as well. 

In addition, the AC phase may be modified as well, as 
the dark matter annihilation rate in this phase is degen- 
erate in the parameter (crw)/mDM- locco et al. fl8i| find 
that the duration of the AC phase may change by almost 
50% if the dark matter mass is changed by a factor of 
2. The effect of different neutralino masses is therefore 
uncertain and should be explored in more detail. We 
will however assume that the general behaviour involv- 
ing adiabatic contraction in the minihalo is still similar, 
such that the calculations below are approximately cor- 
rect also for different neutralino masses. 



VI. COSMIC CONSTRAINTS ON LIGHT DARK 
MATTER 

Obse rvati ons of 511 keV emission in the center of our 
Galaxy |l06l | provi de recent motivation to models of light 
dark matter |l07l |. Such observational signatures can 
be explained assuming dark matter annihilation, while 
other models still have difhculties reproducing the ob- 
servations Q. The model assumes that dark matter an- 
nihilates into electron-positron pairs, which in turn an- 
nihilate into 511 keV photons. Direct annihilation of 
dark matter into gamma-rays or neutrinos is assumed 
to be suppressed to avoid the gamma-ray constraints 



and to ensure a sufficient positron production rate. It is 
known that electron-positron annihilation occurs mainly 
via positronium- fo rmat ion in our galaxy 108]. In addi- 
tion, it was shown '109] that dark matter annihilation to 
electron-positron pairs must be accompanied by a contin- 
uous radiation known as internal bremsstrahlung, arising 
from electromagnetic radiative corrections to the dark 
matter annihilation process. 

Motivated by these results, it was proposed that inter- 
nal bremsstrahlung from dark matter annihilation may 
be responsible for the gamma-ray background at ener- 
gies of 1-20 MeV [98] . Conventional astrophysical sources 
cannot explain the observed gamma-ray background at 
these frequencies [l3| ■ A comparison of the observed and 
predicted background below 511 keV yields constraints 
on the dark matter particle mass [93j . Here we examine 
whether and how this scenario is affected if dark stars 
form in the early universe. We use a thermally aver- 
aged cross section {av) ~ 3 x 10~^^ cm'^ s~^ to account 
for the observed dark matter density [1]. This implies 
that {av) is velocity- independent (S-wave annihilation). 
While Boehm et al. argue that S-wave annihilation 
overpredicts the flux from the galactic center, others ar- 
gue that it is still consistent [H, [sl] . The cross-section 
adopted here is well-within the conservative constraints 
of Mack et al. [9l|. The effect of light dark matter an- 
nihilation on structure formation in the early univers e 
has been studied in various works, e. g. [llOl lllll Ill2| . 
Constraints from upcoming 21 cm observations have been 
explored by Furlanetto et al. JH] and Valdes et al. ^lisj, 
while constraints from background radi ation have been 
considered by Mapelli and Ferrara [ll4| . The effects of 
early dark matter halos on reionization ha ve b een ad- 
dressed recently by Natarajan and Schwarz [ll5j ]. 

As in the previous section, we point out that signifi- 
cant uncertainties are present when considering dark star 
models for different dark matter masses, as this question 
is largely unexplored. In particular, we emphasize that 
no capturing phase will be present for light d ark m atter, 
as shown in the work of Spergel and Press [l05j ]. An- 
other uncertainty is the question whether to adopt self- 
annihilating dark matter (i. e. Majorana particles) or 
particles and antiparticles of dark matter. In the calcu- 
lations below, we assume that light dark matter is self- 
annihilating. Otherwise, our results would be changed 
by a factor of 0.5. 

A. 511 keV emission 

The expected X-ray background from 511 keV emission 
is calculated from Eq. ([5]). The volume emissivity of 511 
keV photons is given as 

Py ^ 5 ((1 + z)v - j/5ii) 511 keVa5ii {av)n^y[C^ii, 

. . (^6) 

where {av) denotes the thermally-averaged annihilation 
cross section, a^n is the fraction producing an electron- 
positron pair per dark matter annihilation process and 
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FIG. 7: The predicted X-ray background due to 511 keV 
emission for different dark matter particle masses. Solid 
lines; Enhanced signal due from adiabatic contraction, dot- 
ted lines: Conventional NFW profiles. The observ ed X -ray 
background from the HEA P ex periments (squares) [ll6| and 
Swift/BATSE (triangles) [TTtI I is shown as well. The com- 
parison yields a lower limit of 10 MeV on the dark matter 
mass for the adiabatically contracted profiles, and 7 MeV for 
standard NFW halo profiles. 



ly^ii the frequency corresponding to 511 keV. In our 
galax y, this process happens via positronium formation 
[lOSj . and we assume that the same is true for other 
galaxies. In 25% of the cases, positronium forms in a sin- 
glet (para) state which decays to two 511 keV photons, 
whereas 75% form in a triplet (ortho) state which decays 
into a continuum. We thus adopt 01511 = 1/4 for 511 keV 
emission. C511 refers to the dark matter clumping factor, 
which is still highly uncertain. The main uncertainty is 
due to the adopted dark matter profile and the assump- 
tions regarding substructure in a halo [ol, [sl, [ol, [o^ . 

Here we use the clumping factor for a NFW dark mat- 
ter profile [97] which has been derived by Ahn and Ko- 
matsu [93i, ^981], as to our knowledge, no calculations of 
dark star formation are available for other dark matter 
profiles. For z < 20, it is given as a power-law as 



10^ — 10""^° Mq, corresponding to their calculation of the 
filtering mass. We also calculate the filtering mass ac- 
cording to the approach of Gnedin and Hui [ll] , Gnedin 
[45'], but obtain somewhat lower masses, with ~ 10^ Mq 
halos at the beginning of reionization and ~ 3 x 10'' Mq 
at the end [l^] . This is also in agreement with numerical 
simulations of Greif et al. [59| which find efficient gas col- 
lapse in halos of 10^ Mq. The discrepancy may also be 
due to their different reionization model, which assumes 
reionization to start at redshift 20. 

We recall that the clumping factor can be considered as 
the product of the mean halo overdensity, the fraction of 
collapsed halos above a critical scale and the mean " halo 
clumping factor" that describes dark matter dumpiness 
within a halo. To take into account that electron-positron 
annihilation occurs only in halos above the filtering mass 
Mf , we thus rescale the results of Ahn and Komatsu 

faalii as 



^ _ /coll(Mf) 

L'511 — 7 777 rL^DM/onh, 

JcolU-lWrninj 



(18) 



where the factor /cnh is given from Eq. p^ . As above, 
we assume that the halo mass must be larger than the 
filtering mass, and lower than the critical mass scale Mc 
that corresponds to virial temperatures of 10^ K. For 
comparison, we will also calculate 511 keV emission with 
/cnh = 1- We note that the resulting background will 
be somewhat lower than the result of Ahn and Komatsu 
(93| . as we adopted asn — 1/4 and include only halos 
above the filtering mass scale in the clumping factor. In 
Fig. [3 we compare the results with the obser ved X -ray 
background from the H EAO -experiments [12^ [ll6j | and 
SWIFT ;i29]/BATSE [HI observations 11%. In the 
standard NFW case, we find a lower limit for the dark 
matter particle mass of 7 MeV. For the case with adia- 
batically contracted profiles due to dark star formation, 
we find a slightly higher lower limit of 10 MeV. This is 
because the enhancement is effective only for frequencies 
hiy < 100 keV, where the observed background is signifi- 
cantly larger than at 511 keV, where Ahn and Komatsu 
[q^ I obtained their upper limit. 



Cdm = CDM(0)(l + z)- 



(17) 



where Cdm(O) describes the clumping factor at redshift 
zero and /? determines the slope. For a NFW profile 
[13, Cdm(O) ~ 10^ and (3 ~ 1.8. The effects of dif- 
ferent clumping factors will be explored in future work 
[118*1 . Ahn and Komatsu (93l . [98j included contributions 
from all halos with masses above a minimal mass scale 
Mniin, which was given as the maximum of the dark 
matter Jeans mass and the free-streaming mass. This 
approach assumes instantaneous annihilation of the cre- 
ated electron -pos itron pairs. As pointed out by Rasera 
and Teyssier |119l | , the assumption of instantaneuos anni- 
hilation is only valid if the dark matter halo hosts enough 
baryons to provide a sufficiently high annihilation proba- 
bility, postulating this to happen in halos with more than 



B. Internal Bremsstrahlung 

The internal bremsstrahlung is calculated according to 
the approach of Ahn and Komatsu [1^ . The background 
intensity is given by Eq. ([S]), with a proper volume cmis- 
sivity 



4a g{iy) 



(19) 



where a ~ 1/137 is the finestructure constant and glv) 
is a dimensionless spectral function, defined as 



(20) 
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FIG. 8: The predicted gamma-ray background due to 
bremsstrahlung emission for different dark matter particle 
masses. Solid lines: Enhanced signal due from adiabatic con- 
traction, dotted lines: Conventional NFW profiles. The lines 
overlap almost identically, as the main contribution comes 
from redshift zero, where the clumping factor is large and 
dark stars are assumed not to form. The observed ga mma - 
ray background from the HEAP experiments (squares) IllSH . 
Swift/BATSE (triarigles) pTj . COMPTEL (crosses) [l20t | 
and SMM (plusses) [121 1] is shown as well. 
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FIG. 9: The upper limit of X-ray radiation due to dark 
star remnants. The observed X-r ay background from the 
HEAD experiments (squares) [HI and Swift/BATSE (tri- 
angles) jll7l | is shown as well. Only for very low dark matter 
particle masses, the upper limit is somewhat higher than the 
observed background. However, the actual contribution may 
be lower by some orders of magnitude (see discussion in the 
text). 



with s = 4mDM('7iDM — hv). As Ahn and Komatsu [qS^ 
pointed out in a 'Note added in proof, bremsstrahlung is 
emitted in all dark matter halos, regardless of the bary- 
onic content. There is thus no need to consider any shift 
in the minimal mass scale, the only thing to take into ac- 
count is the enhancement of annihilation due to the AC 
profiles. The clumping factor Cbroms is thus given as 

Cbrcms = C'DM/cnh, (21) 

where /onh is given by Eq. (|13p . In Fig. [51 we compare the 
results with the observe d ga mma-ray background from 
the HEAO-experiments [116} and SWIFT/BATSE ob- 
servations IutII . as well as SMM [l3l| jl2l| and Comptel 
|132| | data [l20j | . We find that the signal is almost un- 
changed in the model taking into account dark star for- 
mation. The reason is that dark stars form mainly at high 
redshifts, in the range where Mp < M^, while the domi- 
nant contribution to the background comes from redshift 
zero. Our results agree with Ahn and Komatsu [os} . 

C. Emission from dark star remnants 

As in ijV CI we consider a scenario where dark stars 
explode at the end of their lifetime and dark matter an- 
nihilation products in their remnants contribute to the 
cosmic background. For simplicity, we consider 511 keV 
emission only, which is also mostly sensitive to modifi- 
cations at high redshift. Again, we assume that dark 
stars form in halos between the filtering mass Mp and 
the mass corresponding to a virial temperature of 10"* K, 



Mc- In this case, the volume emissivity is given as 



= 5((l + z)i/-i/5ii)511 keVasii /r/a 

„ dfcowiMp) . 

X ^DMjcorc -Jl , (,^^) 



where riDM is the mean proper number density of dark 
matter particles and dfcoii/dt can be evaluated from Eq. 
([5|). The model-dependent factor fr determines which 
fraction of the dark matter in the star will be left in 
the remnant, we adopt = 1 to obtain an upper limit. 
The factor /a determines the fraction of the remaining 
dark matter which actually annihilates, which we set 
to /a = 1 as well. As in Wl A\ a^n = 1/4 is the 
fraction of electron-positron annihilations per one dark 
matter annihilation process, corresponding to annihila- 
tion via positronium formation. In Fig. [9l we compare 
the results with the obser ved X-ray background from the 
HEAO-experiments [ull and SWIFT/BATSE observa- 
tions [117,] . 

For dark matter particle masses below 30 MeV, the 
upper limit found here is higher than the observed back- 
ground. Again, as discussed in § IV C[ there are signif- 
icant uncertainties regarding the question whether this 
high contribution can be reached, both due to uncer- 
tainties in the dark star models, which have not been 
explored for light dark matter, as well as the impact of 
a supernova explosion on the dark matter cusp. These 
possibilities should be addressed further in future work. 
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VII. SUMMARY AND DISCUSSION 

In this work, we have examined whether the suggestion 
of dark star formation in the early universe is consistent 
with currently available observations. We use these ob- 
servations to obtain constraints on dark star models and 
dark matter properties. From considering cosmic rcion- 
ization, we obtain the following results: 

• Dark stars with masses of the order 800 Mq as sug- 
gested by Freese et al. [l^ can only be reconciled 
with observations if somewhat artificial double- 
reionization scenarios are constructed. They con- 
sist of a phase of dark star formation followed by 
a phase of weak Pop. II star formation and a final 
star burst to reionize the universe until redshift 6. 

• The same is true for dark stars in which the number 
of UV photons is significantly increased due to dark 
matter capture, as suggested by locco et al. p^ . 

• It appears more reasonable to require that dark 
stars, if they were common, should have similar 
properties as conventional Pop. Ill stars. For MS- 
dominated models, this requires that typical dark 
star masses are of order 100 Mq or below. For 
CD models it requires a dark matter density above 
10^^ — 10^^ GeV cm^'^ if a spin-dependent elastic 
scattering cross section of i x 10"'^^ cm^ is assumed 

m- 

• Alternatively, it may imply that the elastic scatter- 
ing cross section is smaller than the current upper 
limits, that the dark matter cusp is destroyed by 
mergers or friction with the gas or that the star is 
displaced from the center of the cusp. 

• A further interpretation is that dark stars are very 
rare. This would require some mechanism to pre- 
vent dark star formation in most minihalos. 

• However, if the double-reionization models are ac- 
tually true, it would indicate that dark stars form 
only at redshifts beyond 14, which makes direct ob- 
servations difficult. 

• We also note that 21 cm observations may either 
confirm or rule out double-reionization models. 

We have also examined whether the formation of dark 
stars and the corresponding enhancement of dark matter 
density in dark matter halos due to adiabatic contrac- 
tion may increase the observed X-ray, gamma-ray and 
neutrino background. Here we found the following re- 
sults: 

• For massive dark matter particles, direct annihi- 
lation into gamma-rays provides significant con- 
straints for masses less than 30 GeV. 



• For massive dark matter particles, the contribution 
from direct annihilation into neutrinos is well below 
the observed background. 

• In light dark matter scenarios, the 511 keV emis- 
sion is significantly enhanced below frequencies of 
100 keV in the observers restframe. For a certain 
range of parameters, this emission may even form 
a significant contribution of the total X-ray back- 
ground. In this case, we derive a lower limit of 10 
MeV for the dark matter particle mass (while we 
find 7 MeV for standard NFW profiles). 

• In light dark matter scenarios, the background 
radiation due to internal bremsstrahlung is not 
affected significantly from adiabatic contraction at 
early times, as the main contribution comes from 
low redshift. 

• Both for light and massive dark matter particles, 
the annihilation products in the remnants of dark 
stars may provide significant contributions that 
may be used to constrain such models in more 
detail. However, whether this contribution can 
be reached is highly model-dependent and relevant 
questions regarding the death of dark stars has not 
been explored in the literature. 

Future observations may provide further constraints on 
this exciting suggestion. Small-scale 21 cm observations 
may directly probe the HII regions of the first stars and 
provide a further test of the luminous sources at high red- 
shift, and extremely bright stars might even be observed 
with the James- Webb telescope, if they form sufficiently 
late. With this work, we would like to initiate a discus- 
sion on observational tests and constraints on dark stars, 
which may tighten theoretical dark star models and pro- 
vide a new link between astronomy and particle physics. 
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